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Abstract. We study the existence and stability of periodic travelling-wave solutions for 
generalized Benjamin-Bona-Mahony and Camassa-Holm equations. To prove orbital sta- 
bility, we use the abstract results of Grillakis-Shatah-Strauss and the Floquet theory for 
■ periodic eigenvalue problems. 
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1. Introduction. 

Consider the following equation 



^> ut + (a(u)) x - u xxt = lb' («)y + b(u)u xx \ (1-1) 

m 



where a, b : R — > R are smooth functions and a(0) = 0. In this paper we study 
the problems of the existence and stability of periodic travelling-wave solutions u = 
(p(x — vt) for (1.1). It is easy to see that whatever a, b be, the equation for <p has 
no dissipative terms. Hence, any travelling-wave solution of (1.1) is determined from 
Newton's equation which we will write below in the form ip' 2 = U(ip). Therefore 
by using the well-known properties of the phase portrait of Newton's equation in 
the (99, <y9 / )-plane, one can establish that under fairly broad conditions, (1.1) has at 
least one three-parameter family of periodic solutions ip(y) = ip(v, c\\ ip ; y) where c\ 
is a constant of integration and ipo = mimp (see Proposition 1). The parameters 
v, C\ determine the phase portrait whilst ipo serves to fix the periodic orbit within. 
Moreover, if T = T(v, c 1; <p ) is the minimal (sometimes called fundamental) period of 
ip, then ip has exactly one local minimum and one local maximum in [0, T). Therefore 
(p' has just two zeroes in each semi-open interval of length T. By Floquet theory, 
this means that <p' is either the second or the third eigenfunction of the periodic 
eigenvalue problem obtained from the second variation along ip of an appropriate 
conservative functional M(u). If the first case occurs, then one can use the abstract 
result of Grillakis-Shatah-Strauss (|25j) to prove orbital stability whenever d(v) = 
(d 2 /dv 2 )M(ip) is positive. 

In the periodic case we deal with, it is not always so easy to determine the sign of 
d(v). To overcome this problem, we first establish a general result (see Proposition 



1 



6) expressing d(v) through some special line integrals along the energy level orbit 
{H = h} of the Newtonian function H(X,Y) = Y 2 — U(X) which corresponds to 
ip. When a, b are polynomials, these are complete Abelian integrals and one can 
apply methods from algebraic geometry (Picard-Fuchs equations, etc.) to determine 
the possible values of v, c\ and ip where d(v) changes sign. Let us mention that 
even for v and c\ fixed, the sign of d(v) might depend on the amplitude of (p (ruled 
by tpo) as shown in Proposition 8. In this connection, we calculate explicitly the 
main term of d(v) in the case of arbitrary small-amplitude periodic solutions (p of 
(1.1), see formula (7.7). It is shown that the main term depends on the first two 
isochronous constants related to the center (X , 0) into which the orbit (ip, ip') shrinks 
when e = m&xp — miwp — > 0, and on X itself as well. 

We apply our results to prove orbital stability for several particular examples. 

Theorem I. (The modified BBM equation). Let a(u) = 2uu + [3u 3 , b(u) = 0, (3 > 
and u = tp(x — vt) where v > 0, <p(y) = <p(v, 0; ipo; y) be a periodic travelling-wave 
solution of (1.1) which does not oscillate around zero. Then ip is orbitally stable in 
any of the cases: 

(i) 3v 2 - 8c; 2 > 0; 

(ii) 3v 2 — 8u 2 < 0, 2v 2 — 2uv — u 2 > and the periog of ip is sufficiently large. 

Theorem II. (The perturbed single-power BBM equation). Let a(u) = (3u 2 , b(u) = 
•yfiu, (3 > 0, and let u = p(x — vt) where v > 0, ip(y) = <p(v, 0; po; y) be a periodic 
travelling-wave solution of (1.1). Then ip is orbitally stable for small \^\. 

Theorem III. (The perturbed single-power mBBM equation). Let a(u) = (3u 3 , 
b{u) = 7/3-u 2 , (5 > 0, and let u = ip(x — vt) where v > 0, ip(y) = <p(v,0;ipo;y) be a 
periodic travelling-wave solution of (1.1) which does not oscillate around zero. Then 
<p is orbitally stable for small \ 

Theorem IV. (Small- amplitude waves of the perturbed BBM equation.) Let a(u) = 
2uu + \u 2 , b(u) = 7<?('u) and u = ip(x — vt) where v > 0, ip(y) = ip(v,ci,ipo',y) 
be a periodic travelling-wave solution of (1.1) having a small amplitude. Then ip is 
orbitally stable for small \ j\ and (co/v,ci/v 2 ) taken in appropriate domain £1 C M 2 . 

Theorem V. (Small- amplitude waves of the perturbed mBBM equation). Let a(u) = 
2uju + /3u 3 , b(u) = 7 <?(«), P > and u = <p(x — vt) where v > 0, <p(y) = ip(v, 0; ip ] y) 
be a periodic travelling-wave solution of (1.1) which has a small amplitude and does 
not oscillate around zero. Then ip is orbitally stable for 3v 2 — 8u 2 > and small \ 

We point out that, unlike the other cases, in Theorem IV the constant of integration 
ci is not fixed, therefore we consider the whole family of small-amplitude waves. The 
explicit expression of Q is given in the proof. 

Let us mention that for a(u) = 2ku + |w 2 and b(u) = u, equation (1.1) becomes 
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the well-known Camassa-Holm equation 



ttf ~\~ 2ku x ~\~ 3uU x — t^xxt — ^^x^xx ~t~ ^^xxx' (1*^) 

Equation (1.2) was derived as a bi-Hamiltonian generalization of the Korteweg-de 
Vries equation [23] and later Camassa and Holm [10] recovered it as a water-wave 
model. The Camassa-Holm equation is locally well-posed in H s for s > |. Moreover, 
while some solutions of equation (1.2) are global, others blow up in finite time (in 
both the periodic and non-periodic cases) [TTj [TJl Q21 [EH d51 ESI H2]- The solitary- 
waves of Camassa-Holm equation are smooth in the case k > and peaked for k = 0. 
Their stability is considered in [HI HH UHl [26l [271 [36l EZ] . 

For a{u) = 2ku + ^u 2 and b{u) = ju, equation (1.1) serves as a model equation for 
mechanical vibrations in a compressible elastic rod [2"Tl |2"2"] . Some problems such as 
well-posedness and blowing-up in this case were studied in [HJ S21 SI] , and stability 
of solitary waves was investigated in [191 E3] • 

If b(u) = in (1.1), one obtains the generalized Benjamin-Bona-Mahony (gBBM) 
equation for surface waves in a channel [B] . All solutions are global and their solitary 
waves are stable or unstable depending on a(u) [301 1^0] - 

For equation (1.1), the well-posedness and stability of solitary waves in the case 
a(u) = 2ku + 2 ^m p+1 and b(u) = u p are studied in [26J. For k = 0, equation (1.1) 
admits peaked solitary wave solutions, which are stable (see [27]). 

The existence and stability of periodic travelling waves for nonlinear evolution 
equations has received little attention. Recently Angulo, Bona and Scialom in [1] 
developed a complete theory of the stability of cnoidal waves for KdV equation. The 
solution u(x, t) = ip c (x — ct) of KdV satisfies the equation 

fl + ifl ~ CL Pc = A Vc , 
where A iPc is an integration constant. An explicit form for (p c in the periodic case is 



Vc(0 =P2 + (ft - &)cn 2 (y^y/k; 



where cn is the Jacobi elliptic function and the following relations take place: 
ft < /3 2 < ft, k 2 = ft + ft + ft = 3c, A Vc = -iE^-ftft, 



Since cn(u + 2K) = —cn(u) where K = K(k) = Jq — ^ == is the complete 

elliptic integral of the first kind, then (p c has the fundamental (i.e. minimal) period 
T Vc given by 

Moreover, T [fic depends on the speed c and satisfies the inequality 

(2vr) 2 



T > 



c 2 + 2A 
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Using the abstract results of Grillakis, Shatah and Strauss (adapted to the periodic 
context), the authors proved in [I] that the cnoidal waves with mean value zero are 
orbitally stable. 

Other new explicit formulas for the periodic travelling waves based on the Jacobi 
elliptic function of type dnoidal, together with their stability, have been obtained by 
Angulo jH |3] for the nonlinear Schrodinger (NLS) equation iu t + u xx + \u\ 2 u = 0, 
modified KdV equation u t + 3u 2 u x + u xxx = and Hirota-Satsuma system 

u t — a(u xxx + 6uu x ) = 2bvv x 
v t + v xxx + 3uv x = 0. 

For NLS the solutions are of the form u(x,t) = e lwt p w (x) where <p w {x) is a periodic 
smooth function with period L > 0. The solution (p w is of the form 

where 771 and the modulus k depend smoothly on w. The orbit fl<p w , 

n Vw = {e ie ip w (- + y), (y,0) 6lx [0,2tt)} 

generated by the dnoidal wave ip w is stable by perturbation of periodic function with 
period L and nonlinearly unstable by perturbation of periodic function with period 
2L. In all these works it was necessary to use an elaborated spectral theory for the 
periodic eigenvalue problem, 

^ + [p-n{n + l)k 2 sn 2 {x;k)]V = 0, 
*(0) = V(2K(k)), *'(0) = *'(2A'(ife)), ' 

with specific values of n G N. We will also make use of systems similar to (1.3). 

The present paper is organized as follows. In section 2 we formulate and sketch 
the proof of a local well-posedness result for the equation (1.1) in periodic H s spaces. 
In Section 3 we prove the existence of periodic travelling waves of a given (admis- 
sible) period and study their properties. In Section 4 we prove the orbital stability 
result under some hypotheses (see Assumption 1). In Section 5 several particular 
examples are considered. For most of them, we determine ip explicitly and show that 
if' is the second eigenfunction of the respective periodic eigenvalue problem. We also 
determine the sign of d(v) to outline the cases satisfying Assumption 1. In Section 
6 perturbation theory is applied to prove the orbital stability of periodic travelling 
waves for generalized Benjamin- Bona- Mahony and generalized Camassa-Holm equa- 
tions (the case of small b in the right-hand side of (1.1)). In Section 7 we study the 
small-amplitude periodic travelling waves of (1.1) and determine the sign of d(v) for 
them. 

2. Local well-posedness. 
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In this section, we discuss the local well-posedness of the Cauchy problem for equation 
(1.1). We begin by introducing some notation and by recalling related definitions we 
shall use throughout the paper. 

Let V = C™ r denote the collection of all functions which are C°° and periodic 
with a period T > 0. The topological dual of V will be denoted by V. If * G V 
then we denote by ^f(f) = (\&, f) the value of ^ at /. Define the functions Qk{x) = 
exp(2nikx/T), k G Z. The Fourier transform of ^ G V is the function \P : Z — > C 
defined by ty(k) = ©-&)• If ^ is a periodic function with a period T, we have 

l r T 

tf(Jfe) = — / * (x) exp(-2nikx/T)dx. 
T Jo 

For s £l, the Sobolev space H s ([0,T]) is the set of all / G V such that 



11/112 = T E (l + |A:| 2 ) s |/(A;)| 2 <oo. 

k=— oo 

Certainly, H s ([0,T]) is a Hilbert space with respect to the inner product 



(f,9)s = T E (l + \k\ 2 ) s f(k)g(k). 

k=— oo 

Since if s ([0,T]) C L 2 ([0,T]) for every s > 0, we obtain via Plancherel identity that 
for every n G N, 

n 

II/IG = Ell/ a) ll 2 - 

J"=0 

where represents the jih derivative of / taken in the sense of V'. Moreover, 
Sobolev's lemma states that if s > I + |, then H 2 ([0,T]) <^-> Cp er where 

Cper = {/ G : f ij) is periodic with a period T for j = 0, I}. 

One can prove the following result about equation (1.1). 

Theorem 1. Assume that a,b G C m+3 (R), m > 2. Given uq G -£P, | < s < m, there 
exists a maximal t > and a unique solution u(x,t) to (1.1) such that 

ueC([0,t o ),H s )nC 1 ([0,t o ),H s - 1 ). 

Moreover, the solution depends continuously on the initial data. 
Proof. Take u G H s and let 

A{u) = b(u)d x , f(u) = (1 - d 2 x )- l [b(u)u x - d x {\b'{u)u 2 x + a(u))}. 
Using the above notations, one can rewrite equation (1.1) in the following form: 

u t + A(u)u = f(u). 
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In a similar way as in Theorem 2.2 in [27] (dealing with the non-periodic case), we 
have 

(1) A(u) is quasi-m-accretive, uniformly on the bounded sets in H s ~~ l . Moreover, 
A(u) G L(H S , iP -1 ) (where L(X, Y) is the space of all linear bounded operators from 
X to Y, L(X) = l[x,X)) and 

||(A(«i) - A(u 2 ))u 3 \\ s -1 < (JLl\\ui - 'U2||s-li|«3||a- 

(2) Define A = (1 - d 2 ) 1 ' 2 , B(u) = [A,6(m)9JA~ 1 for u G H s , where [A, M] denotes 
the commutator of A and M. Then B(u) G L^H 8 ^ 1 ) and 

\\{B{u x ) - B(tt 2 ))u 3 || s _i < /i 2 ||wi -UzlUlualla-i, Ui,"2 G m 3 G 

(3) /(«) is bounded on the bounded sets in H s and satisfies 

||/(«i) - f{u 2 )\\ s < //allwi - W2IU «i,«2 e 

||/(tti) - /(u 2 )|| a _i < m\\u\ — «i, 1*2 e # s-1 . 

Applying Kato's theory for abstract quasilinear evolution equations [22], we obtain 
the local well-posedness of the equation (1.1) in H s , for | < s < m. The solution 
u(x,t) belongs to C([0,t ), H s ) D C\[0,t ), H 3 - 1 ). 

3. Periodic travelling-wave solutions. 

We are looking for a travelling-wave solution of (1.1) of the form u(x, t) = tp(x — vt). 
We assume that ip is smooth and bounded in R. The following two cases appear: 

(i) ip' 7^ in R and ip- < ip < ip + (corresponding to kink- wave solution); 

(ii) — for some Denote ipo = <p(£), <p>2 = ¥>"(£)• 

Below we will deal with the second case. Replacing in (1.1) we get 

-vif' + (a(ip))' + vif'" = (b'{y)^- + b{ip)A . (3.1) 
By integrating (3.1) twice, one obtains 

-vip + a(<p) + vif" = b'(ip)^- + b((p)<p" + ci, (3.2) 

ip' 2 v ff 

[ v -b(<p)} — = c 2 + a<p + -ip 2 - A(<p), A(<p)= a(s)ds, (3.3) 
II Jo 

with some constants ci,c 2 . In case (ii), one has respectively 
Ci = a(ipo) - vip + (v- b(<p ))(p 2 , 

c 2 = A(ip ) + \vifl - <p Q a(jp ) - (v - b(jpo))<p <p 2 = A(ip ) - \vif% - ci<p . 
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Next we are going to look for periodic travelling-wave solutions ip. Consider in the 
plane (X, Y) = (ip, ip') the Hamiltonian system 



X = [v-b(X)]Y = H Y , (3 n 

Y = \b'{X)Y 2 - a{X) + vX + Cl = -H x 



with a Hamiltonian function 

H(X, Y) — [v — b(X)}^ + A(X) - V -X 2 - Cl X. 

Then (3.3) becomes H(ip, ip') = c<i and the curve s — > (ip(s-so), <p'(s — so)) determined 
by (3.3) lies on the energy level H = C2 of the Hamiltonian H(X,Y). Within the 
analytical class, system (3.4) has periodic solutions if and only if it has a center. 
Each center is surrounded by a continuous band of periodic trajectories (called period 
annulus) which terminates at a certain separatrix contour on the Poincare sphere. As 
far as the straight lines X = X where v — b(X ) = are (unions of) trajectories, the 
critical points of center type of (3.4) are given by the critical points on Y = having 
a negative Hessian. These are the points (A , 0) where: 

ci + vX - a(X ) = 0, [v- b(X )\ [v - a\X )] < 0. (3.5) 

(For simplicity, we will not consider here the case of a degenerate center when the 
Hessian becomes zero.) 

The above considerations lead us to the following statement. 

Proposition 1. Let c\ and v be constants such that conditions (3.5) are satisfied for 
some Xq G K. Then there is an open interval A containing Xq such that: 

(i) For any ip G A, ip < X , the solution of (1.2) satisfying 

<AU = <Po, V (U =0, <p (0 = — r , 

v - b(ipo) 

is periodic. 

(ii) If ipi G A ; ipi > X is the nearest to X solution of H(X, 0) = H(ip ,0), then 
Po<(p< <P>\- 

(iii) If T is the minimal period of ip, then in each interval [s, s + T) the function 
ip has just one minimum and one maximum (equal to ip and <p>\, respectively) and is 
strictly monotone elsewhere. 

Proof. Part (i) follows from the analysis already done above. Parts (ii) and (iii) 
follow from the fact that H(X, Y) is symmetric with respect to Y and because for 
any fixed X the equation H(X, Y) = const has just two solutions Y, —Y. 

Example. Assume that a, b are polynomials and dega = 2, degb < 1. For ci, 
v properly chosen, the quadratic equation in (3.5) will have two distinct real roots 
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X\ < X 2 . It is easy to see that if v — b(X) ^ in [Xl,X 2 ], then (3.5) holds for just 
one of these roots. If a"V < and v — b(X) vanishes at X 1; then (3.5) holds for X 2 , 
and vice versa. Finally, if a"b' < and v — b(X) has a root within (Xi,X 2 ), then both 
X\ and X 2 satisfy (3.5). To summarize: the phase portraits of all quadratic systems 
(3.4) (including the exceptional case dega < 1, deg 6=1) satisfying (3.5), that is 
having a center, are divided into 7 topologically different classes, see e.g. 



Remark 1. Below, we shall denote A~ = A n {(-oo,X )}, A+ = A n {(X ,oo)}. 
It is possible that (3.3) also has periodic solutions for initial values (po far from Xq. 
To study them, it is needed to specify the functions a and b in (1.1). 

Let us denote 

2c 2 + 2cis + vs 2 - 2A(s) 2A(ip ) - vipft - 2ciy? + 2cis + vs 2 - 2A(s) 



v — b(s) v — b(s) 



U(s) 



Then for (p < (p < (fx one can rewrite (3.3) as (p'(cr) = JU(cp(a)). Integrating the 
equation along the interval [£, s] C [£,, £ + T/2] yields an implicit formula for the value 
of <p(s): 

r<p(s) Arc 

\ -TT^ = S -^ se[Z,d + T/2}. (3.6) 

For s G [£ + T/2, £ + T] one has ip(s) = <p(T + 2^ — s). We recall that the period 
function T of a Hamiltonian flow generated by H = ^Y 2 — ^U(X) = is determined 
from 

f T , r dX m dX , . 

T= dt=f -^ = 2 -== (3.7) 
Jo Jh =o Y Jipo Ju(X) 

This is in fact the derivative (with respect to the energy level) of the area surrounded 
by the periodic trajectory through the point (<p , 0) in the (X, Y) = (ip, y/)-plane. 

Consider the continuous family of periodical travelling- wave solutions {u = (p(x—vt)} 
of (1.1) and (3.3) going through the points (<p, <p') = (<po,0) where ip E A - . For any 
<p G A - , denote by T = T(ip ) the corresponding period. One can see (e.g. by using 
formula (3.7) above) that the period function <p — > T(ip ) is smooth. To check this, 
it suffices to perform a change of the variable 

X = ^Z_fl s + ^±^ (3 .8) 

in the integral (3.7) and use that 

£% ) = = 0. (3.9) 

Also, it is not difficult to verify (see Section 7) that 



T((p ) -* T = 2tt a 



\ a'(X )-v 



b(X ) 

as ip t X 
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Conversely, taking v, c\ to satisfy the conditions of Proposition 1 and fixing T in a 
proper interval, one can determine ipo and p\ as smooth functions of v, C\ so that the 
periodic solution p given by (3.6) will have a period T. The condition for this is the 
monotonicity of the period. 

Definition 1. We say that the period T = T(ip ) is not critical provided that 

If the period T(p ) is not critical for any p G A_, then the period function is strictly 
monotone along the period annulus and its range is an open interval / having T as 
an endpoint. 

Proposition 2. Let A be a period annulus of (3.4) which surrounds a nondegenerate 
center and has a monotone period function. Then for any T G I there is a unique 
p G A_ satisfying T(p ) = T. Its derivative ipo with respect to v is determined from 



Mci+W0 -a {v ^i * = U (X (3.10) 
dh JH=h y da jh=k 2y 



l 2 _ 

where 



h = A(p Q ) - -pi - ap . (3.11) 

Proof. We use the implicit function theorem (IFT) and the Gelfand-Leray form (see 
[6], Chapter 3). Let I be the range of the period function along A. For Tel and 
p G A_ denote 

G(v, ci, p ) = G(v, c h h) = T(p ) —T=<f — — T. 

JH=h y 



Then 



/-i / rpi ( \ 9G dh dG d I dx 

opo dpo Oh dh JH=h y 



and the IFT works. Hence 

=, dG • dG , _ x 

From H(x,y) = h we obtain the covariant derivatives 

(v - b (x))y% = !> \(y 2 - + (?- h ^))yy = o. (3.13) 

Therefore, by using Gelfand-Leray form to calculate the derivatives, we get 

^, ,n d f , ,, u , m dG d f dx , 
G{v,c 1: h) = — j (v-b(x))ydx-T, — = — j — , 3.14 
dh Ju=h oh dh JH=h y 

dG d r , . ,, , d r (x 2 + y 2 )dx 

ftT = dh U v + (v ~ b{xMix = Th L h ar- (X15) 
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Hence, by (3.12), (3.14) and (3.15) we obtain 



h±{ - + &+/ )dX =°- (3-16) 

ah JH=h y ah JH=h 2y 

Finally, from (3.11) one obtains h = — ^V^o + [ a (Vo) ~~ V( Po ~ c i]0o- Together with 
(3.16), this implies (3.10). □ 

Remark 2. Obviously, one can formulate a local variant of Proposition 2 concerning 
a given noncritical period T(ipo) only. As far as T'(tpo) ^ 0, the same proof clearly 
goes and no restrictions concerning the period annulus A are needed. 

The perturbation result we establish below will be needed in Section 6. Instead of 
(1.1), consider now a small perturbation of the generalized BBM equation 

u t + (a(u)) x -u xxt = ^\b\u)^ + b(u)u xx ) , ItI « 1 (1-1 7 ) 



X 



and let {^(x — vt)} be the family of corresponding periodic travelling-wave solutions 
going through points (</?q, 0) in the (ip, <//)-plane. Denote the related periods by T((po). 

Proposition 3. Assume that the period T = T(tpo) related to the solution tp°(x — vt) 
of (1.1 ) is not critical. Then: 

(i) There is a smooth function 7 — > (fo(l) defined for small \j\ and satisfying 
<Po(0) = y?o; su °h that the travelling wave solution ^(x — vt) of (1.1 7 ) going through 
the point (vo(t),0) ^ periodical and has a (minimal) period T. 

(ii) max [0i T] |</? 7 - </?°| = 0{l)- 

Proof, (i). The proof relies on the implicit function theorem. Consider system 
(3.4) with b replaced by 76. Since for ip e A - the point (ipo, 0) is not critical for 
(3.4), IFT yields that there is a smooth function </?i = <£>i(ty?o, 7) determined from 
H(ipi,0) = H(ifi O ,0), which takes values in A + and ip < p 1 < ipi. Introducing a 
new variable (3.8) in (3.7) and making use of (3.9), we can rewrite (3.7) in the form 
(7(7, (po) = where 

f/ " S - T, 



G(i, V o) = 2 ( 



with f/i a smooth positive function. We further have 



GM)=0, 

d(f 



7=0 



therefore (i) follows, again by the IFT. 

(ii) Consider system (3.4) with b replaced by 76 and initial data X(£) — v?o(t), 
Y(£) = 0. By the uniqueness and the smooth dependence theorems, the solution y? 7 
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is smooth and therefore a uniformly Lipschitz continuous function with respect to 7 
for I7I small, on each interval [£, £ + T]. □ 

In support to the hypothesis of Proposition 3, we state the following (in fact, known) 
result about the period T = T((p ). 

Proposition 4. (i) There is no critical period provided that a(u) in (1, 1 7 ) is a 
polynomial of degree 2. (ii). There is at most one critical period, ifa(u) is a polynomial 
of degree 3. 

Proof. Consider the period function related to the Hamiltonian system (3.4) (with 
b(X) = 0) as a function T(h) of the energy level H = h. It is well known that 
T'(h) ^ if deg a = 2 [2U] and that T'(h) has at most one zero if deg a = 3 [24J. As 

[a(<p ) - vtpo - ci}- d 



dip dh 
and the coefficient is not zero for (po G A~ according to (3.5), the claim follows. □ 

4. Conservation laws and conditional stability. 

We now turn to our main stability problem. Take v G R and denote by (p v the periodic 
travelling- wave solution u = <p(x — vt) of (1.1). 

Definition 2. The periodic travelling-wave solution ip v of (1.1) is said to be stable, 
if for every e > there is 5 > such that if u G C([0, t ), H 1 ) is a solution of (1.1) 
satisfying ||w(0) — y?„||i < 5, then 



inf ||ti(-,t) — <f(- — r)\\\ < e for every t G [0,t 



0; 



rell 



For e > 0, consider a neighborhood U £ in the set of all translations of (p v defined by 
U £ = {ue H\[Q, T\) : inf \\u - tp v (- - r)||i < e}. 



rell 



Lemma 1. There exists e > and a C 1 -map a :U £ — > K./T such that for all u eW £ , 
one holds (u(- + a(u)), ip' v ) = 0. 

Proof. Standard (for details, see Lemma 4.1 in [U]). 

As usual, an important role in our construction will be played by some functionals 
which are invariant with respect to time t. Equation (1.1) possesses the following 
conservation laws: 



E(u) 



T 







A(u) + -\±u 2 x 



dx, 



Q{u) = 7: I (u 2 + ul)dx, 
I Jo 

V{u) = I udx 
Jo 
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where T is the minimal period of the solution u(x, t) and A'(u) = a(u). Let us denote 
for short M = E + vQ. In terms of E, Q, and V equation (3.2) with ip = ip v reads 

M\<p v ) + Cl = E'(<p v ) + vQ'(<p v ) + Cl V'((p v ) = 0. (4.1) 

Let 

d(v) = M(<p v ). (4.2) 
Then differentiating (4.2) with respect to v we obtain 

d\v) = Q((p v ), 

d{v) = (<?(<p v ), <p v ) = ±Q + V ' v 2 )dx^ . (4 ' 3) 

Consider in L 2 [0,T] the operator H v defined by the formal differential expression 

H v = M"(tp v ) = (6(^)-u)^+6'(^)^9 x +7;-a , (^) + i6' , (^)^ + ^(^)^. (4.4) 

As 7i v (p' v = 0, zero is in the spectrum of H v . We make the following assumption 
concerning H v and d(v): 

Assumption 1. (i) The operator 7i v has a unique negative eigenvalue, a simple 
eigenvalue and the rest of its spectrum is positive, (ii) d(v) > 0. 

Lemma 2. If Assumption 1 holds and y satisfies (Q'(ip v ),y) = (<p' v ,y) = 0, then 
(H v y,y)>0. 

Proof. Differentiating (4.1) with respect to v yields (since c\ does not depend on v) 

H V V = -Q'(<p v ) (4.5) 

and from (4.3) we obtain 

(H v <p v ,<p v ) = -d(v)<0. (4.6) 

Putting y = a±x + Pi, Pi G P, where x is an eigenfunction of Ti v corresponding to 
the negative eigenvalue — Aq and P is the positive subspace of H v , we obtain 

(H v y,y) = -a\\l + (H v pi,pi). 

Write ip v = a x + bo<-p' v + Po, Po G P- From (4.6) we have 

> (Hv<p v ,<p v ) = -a>l\ 2 + (H v p ,Po) 



and 



Hence 



0= -(Q'((p v ),y) = (H v <p v ,y) 
= (-ooAqX + HvPo,aiX+Pi) = -a aiAg + (H v p ,Pi). 

Inj \ 2 \ 2 , /i_y \ \ 2\2 , (^vP0,Pl) 2 

(H v y,y) =-a 1 X + {H v pi,pi)>-a 1 X + -prr r 

,2\2 , ti u, 1 /\ 2 



~ a i A o + T9~9 77T - a i K i- 
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The rest of the proof is the same as in [29J, pages 310-311. 

Proposition 5. If Assumption 1 holds, then for any positive e there exists K > 
such that for any solution of (1.1) satisfying u ElA e , Q(u) = Q(tp v ), one holds 

E{u) - E(tp v ) + Cl V(u) - c x V( Vv ) > K\\u{. + a{u)) - <p v \\\. (4.7) 

Proof. Denote ip — u(- + a(u)) — (p v — fi(ip v — (p") + y where ((p v — <p", y) = 0. By 
the translation invariant property of Q, we have 

Q(<p v ) = Q{u) = Q(<p v ) + (Q'M,^) + 0(||V||?) 

= Q(if v ) + (<p v - <plfi(<p v - rt) +y) + 0(||V||?) 

= Q{<p v )+v\\<p v -rt\\h+om\i)- 

This implies \i = 0(\\ip\\f). Since Q(u) = Q((p v ), the identity M"(<p v ) = TC V and the 
Taylor expansion give 

E{u) - E(<p v ) + ci / il>dx=-(H v y,y) + o(\\ip\\l). 
Jo 2 

From Lemma 1 we have 

= (u(- + a(u)), (f' v ) = (fi{(f v - iff v ) +y + f v , f'v) = (y, <Pv)- 
From the above equality and Lemma 2 we obtain 

E(u) - E{cp v ) + Cl f ^dx > K\\y\\l + o{\\iP\\l). 
Jo 

This estimate together with 

IMIi = - K<p - > IIV'lli- l^lllv-^lli 

yield 

E{u) - E(tp v ) + ci f ipdx > K\\if}\\l 
Jo 

for sufficiently small. This completes the proof. 

Theorem 2. If Assumption 1 holds, then the travelling-wave solution (p v is stable. 

Proof. Suppose that ip v is unstable. Then there exists a sequence of initial data 
u n (Q) G H 1 and n > such that 

\\u n (0) - <f v \\x -> 

but 

sup inf \\u n (-,t)- (p v (- - r) \\ x > r), (4.8) 

t£[0,t }r&R 
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where u n £ C([0,t ); H 1 ) is a solution of (1.1) with initial data u n (0). Let t n £ [0, t ) 
be the first time so that 

inf ||u„(-,f n ) - <Pv(- ~r)\\i = V- 

r&K 

We have 

%(.,g) = %(o))^%„), 

Q(u n (-,t n )) = Q(u n (0)) -> Q(<p„), 

%(-,g) = %(o))^%). 

Choose a sequence ip n £ H 1 such that Q(ip n ) = Q{fv) and ||*0 n — u n {-, t n )\\i — > 0. By 
continuity of £ and V, E(ip n ) — > £ , (</? t) ) and V(VVi) — > V(v?„). From (4.7) we have 

E(ip n ) - E(<p v ) + ci(y(^ n ) - V^)) > X||Vn(- + a(^„)) - <A,||i- 
Therefore | — ip v (• — «(■?/>«)) ||i — > 0, which implies 

\\u n ~ <fv(- - a(ll>n))\\l -> 0. 

This however contradicts (4.8). The proof of Theorem 2 is complete. 

In order to apply Theorem 2 we have, according to Assumption 1, to determine the 
sign of the derivative d{y). Using the same technique as in the proof of Proposition 2, 
we obtain the following expression of d(v) through line integrals. When a(u) and b(u) 
are polynomials (or even rational functions), these are complete Abelian integrals. A 
lot of methods have been developed to investigate Abelian integrals, which could be 
applied here to study the sign of d(v). 

Proposition 6. Assume that the period T = T({p ) = T(h) is not critical where h is 
given by (3.11). Then 

<«> 

with 

W(h) = (— <£ —)(—{ ( x4 + 2x2 y 2 -& 4 ) d A (d r (x 2 + y 2 )dx \ 2 
\dh JH=h y J\dh JH=h y J \dh Jn=h y J 

Proof. We apply again (3.13) and use the Gelfand-Leray form (in both directions) 
to calculate the needed derivatives. As 

■u n 1 / (x 2 + y 2 )dx If , 1 d r , 

d[v) — - <b = - <z> ydx H — <t> x (v — b(x))ydx 

V ; 2J H =h y 2J H =h 2dhJ H =h V V JJy 

one obtains 

77 x h d r (x 2 + y 2 )dx 1 d f r , d r 2 . r/ ^ , \ 

d(v) = — — <p 1 — f ydx + — <t x (v — b{x))ydx 

KJ 2dhJH=h y 2dv\JH=h y dhJ H =h v v ,w J 

hd r (x 2 + y 2 )dx^r (x 2 — y 2 )dx d r (x 2 y x 2 [x 2 — y 2 ) \ 



2 dh JH=h y JH=h 4(t> — b(x))y dh Jn=h \ 2 Ay 

hr (x 2 + y 2 )dx^ld r (x 4 + 2x 2 y 2 - \y 4 )dx 
2 JH=h y 4dh Ju=h y 
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Replacing the value of h from (3.16), we come to the needed formula. □ 
5. Examples. 

Example 5.1: the BBM equation. Consider the BBM equation 

u t + 2uju x + 3uu x — u xxt = 0, w£l (BBM) 

which is a particular case of (1.1) with a(u) = 2uou + \u 2 and 6 = 0. To apply 
Theorem 2 to BBM, we have to verify Assumption 1. Namely, to establish that the 
corresponding operator Ti v has the needed spectral properties (i) and to prove the 
convexity of d(v). 

Let us first mention that (3.5) reduces now to D = (v — 2uj) 2 + 6c x > and: 
V-2U + D 1 / 2 (V-2U-D 1 / 2 v -2u + 2D 1 / 2 \ r 

*.= , A=( , j ,f „>0, 

X = V -\- D " 2 , A^("- 2 "- 2Z)1/2 ,"-% +Z)1/2 ) if „<0. 
By the definition of Ci, C2 and U(s) one obtains in the considered case 

U (s) = ^(v? - s) [s 2 + (ip + 2u - v)s - (2y?o + (2u - v)ip + 2vLp 2 )} 
(s - ip )(ipi - s)(s + </?i + y? + 2uj - v) 



We note that the last equality is a consequence of Proposition 1 (ii) which implies 
that U(ipi) = U((po) = 0. To obtain an explicit formula for the travelling wave ip v , 
we substitute a = y?o + (y?i — <£o)z 2 , z > in order to express the integral in (3.6) as 
an elliptic integral of the first kind in a Legendre form. If v < 0, one obtains 

dz 

/ , = a(s - £) 

JO J(\ _ 7 2W1 _ i-2 7 2\ 



f 1 -z 2 )(l - A; 2 z 2 ) 

where 



V </?i - </?o V?i + 2^o + 2u; - w ' V 4-u 

Thus we get the expression 

<Pv(s) = y?o + (y?i - y?o)sn 2 (a(s - 0; k). (5.1) 
Similarly, in the case v > we obtain (with the same Z) 



^0 ^(1-^)^/2 + ^2) 
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where 



2 V^l-V^O 12,1/2 1 _ hfO + 2( Pl + 2uJ ~ V 



k 2 = Z^21 ^ + A/ 2 = 1, a 



(^o + 2(^i + 2u; — f V 4f 

and the expression for 

<^ v (s) =f + (ipi - ^o)cn 2 (a(s - £); k). (5.2) 

To calculate the period of ip v , we use (3.7) and the same procedure as above. In 
this way we get in both cases 

vi da 2 /-i dz 2K(k) 



T = 2 - 

'vo ^U{a) & J° J(l -z 2 )(l - k 2 z 2 ) « 

We return to the operator T~C V defined by (4.4) which now has the form 

H v = —vdl + v — 2u — 3y?„ (5.3) 

where (p v is determined by (5.1) or (5.2). Take v > and consider the spectral 
problem 

H v ib = \ip, 

V(0)=V(T), f(0) = ij'(T). 

We will denote the operator just defined again by TC V . It is a self-adjoint operator 
acting in H 2 ([0,T}). From the Floquet theory applied to (5.4) it follows [39J that its 
spectrum is purely discrete, 

Ao < Ai < A 2 < A 3 < A 4 < . . . (5.5) 



where Ao is always a simple eigenvalue. If ip n (x) is the eigenfunction corresponding 
to A n , then 

ipo has no zeroes in [0,T]; 

(5.6) 

ip2n+i, ip2n+2 have each just 2n + 2 zeroes in [0, T). 

Proposition 7. The linear operator 7i v defined by (5.3) — (5.4) has the following 
spectral properties for v > 0: 

(i) The first three eigenvalues of7i v are simple. 

(ii) The second eigenvalue ofTL v is Ai = 0. 

Proof. By (3.1), TC v <f' v = 0, hence ip = <p' v is an eigenfunction corresponding to zero 
eigenvalue. By Proposition 1 (iii) ip' has just two zeroes in [0,T) and therefore by 
(5.6) either = Ai < A 2 or Ai < A 2 = or Ai = A 2 = 0. We are going to verify that 
only the first possibility = Ai < A 2 can occur. From the definition of k and a one 
obtains that 

ip + 2<fi + 2u — v = 4va 2 , <fi — <fo — ^vk 2 a 2 . 
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Then using (5.2) we get 

Ti v = —vd 2 + v — 2uj — 3<y9o — 3(y>i — ipo)cn 2 (ax; k) 
= —vd 2 + v — 2uj — 3<y2i + 3((^i — ip )sn 2 (ax; k) 
= -v&l - va 2 [4k 2 + 4 - 12k 2 sn 2 (ax; k)] 
= va 2 [-d 2 -4k 2 -4 + I2k 2 sn 2 {y; k)] = va 2 A 
where y = ax. The operator A is related to Hill's equation with Lame potential 

d 2 

Aw = — —w + {I2k 2 sn 2 (y; k) - 4k 2 - 4]w = 
ay 2 

and its spectral properties in the interval [0, 2K (k)] are well known [21EIE2]- The 
first three (simple) eigenvalues and corresponding periodic eigenfunctions of A are 

Ho = k 2 - 2 - 2y/l - P + 4A; 4 < 0, 

tf; (y) = dn(y; k)[l - (1 + 2k 2 - y/1 - k 2 + 4k A )sn 2 (y- k)} > 0, 
Hi = 0, 

ih(v) = dn (v, k)sn(y; k)cn(y; k) = \{d/ 'dy)sn 2 (y; k), 
H 2 = k 2 - 2 + 2y/l -k 2 + 4k 4 > 0, 

ij; 2 (y) = dn(y; k)[l - (1 + 2k 2 + y/1 - k 2 + 4k i )sn 2 (y; k)]. 

As the eigenvalues of TC V and A are related by A n = va 2 fi n we conclude that for 
v > 0, the first three eigenvalues of (5.3)-(5.4) are simple and moreover Ao < 0, 
Ai = 0, A2 > 0. The corresponding eigenfunctions are ipo(ax), ipi(ax) = const.ip' v {x) 
and ip2 {ax). □ 

What about the sign of d(v), it is easily seen it is positive when u = (see the 
end of this section). For u 7^ 0, the proof is much more long and complicated. It 
can be done by following the procedure we propose in the next subsection. For this 
reason, we are not going to investigate the case to 7^ in the present paper. 

Example 5.2: the modified BBM equation. Our second example is concerned 
with the periodic travelling-wave solutions of the equation 

u t + 2wu x + (3(u 3 ) x - u xxt = (mBBM) 

where w,/? e R and (3 7^ 0. For this case, a{u) = 2wu + f3v? and b(u) = in (1.1). 
For definiteness, we take tp v — tp(x — vt) where v > 0. The Hamiltonian flow in (3.4) 
is generated by 

H(X, Y) = V -Y 2 + I X 4 + ( w - I) X 2 - Cl X. 
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Our plan is to study here the "symmetric" case C\ = 0. The general case could then 
be considered as a perturbation of the symmetric one (at least for c\ small). It is well 
known that there are three cases related to the symmetric Hamiltonian 



Y 2 + ^X 4 
2v 



2 -^-l)X 2 

v 



v>0, 



see e.g. [31] : 

(i) global center: (3 > 0, 2u > v; 

(ii) truncated pendulum: (3 < 0, 2oj > v ; 

(hi) Duffing oscillator: (3 > 0, 2u < v, 

with three topologically different phase portraits. There is one continuous family of 
periodic orbits in cases (i) and (ii) and three families (left, right, outer) in (iii) (we 
advice the reader to draw a picture). This also holds true for all c\ in case (i) and 
for c\ < -^f3~ l {y — 2u) 3 in (ii), (iii). In the symmetric case c\ = we deal with, the 
periodic solutions exist if and only if ip G A - (see Proposition 1) where: 

A" 

A- 

A" 

A- 



-oo, 0) in case (i 



(^) 1/2 ,0) in case (ii); 



- P > ' 



[^) 1/2 ) in case (iii), left family; 



0. 



' - p^ ) 1 / 2 ) in case (iii), right family; 

-oo,-(^) 1/2 ) 
The respective functions U(s) take the form: 



in case (iii), outer family. 



" A p 

U(s) = -lM-s 2 )( 2 ^ - 



p 



Vl - s 2 ) 



U(s) = l(s 2 - V 2 )( V 2 -s 2 ) 
U{s) = Us 2 -rt){rt-s 2 ) 



in case (i) and case (iii), outer family, 
in case (ii), 

in case (iii), left family, 
in case (iii), right family. 



In a similar way as we have done in the (BBM) case, we can use formula (3.6) to 
calculate explicitly ip v . By an appropriate change of the variables in (3.6), one can 
express the solution through standard elliptic integrals (we omit the details). Up to 
a translation of the argument, we have: 







= ip cn(as; 


k), 


a = 


<fv 




= ip sn(as; 


k), 


a = 


fv 




= ip dn(as 


k), 


a = 


fv 




= ipidn{as 


k), 


a = 



2u—v+P(Pq 



2v ' 

! W 
2v ' 

rp~vf 

2v ' 



k 
k 
k 
k 



1 W 






1 w 


2v- 




Uu>- 


-2v+2Ptp% 


w 


Ua>- 


-2v+2P<fl 
a. .'1 



(i) and (iii) outer, 



(iii) left, 
(iii) right. 
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These formulas of a and fe yield the following expressions and range I for the period: 



T = y^vT+FtfM, fee (0,1), /=(2tt 



2td-i; ' 



OO 



T = 2.-i-V^Fif(fc), fee (0,1), /=( 2 7r./^ i -,oo), (hi) left and right 



T = 4 v ^v / 2fe 2 ~^l^(fe), fee(^,l), J=(0,oo), Router. 

Indeed, the formulas just derived imply that T = T(k) is strictly decreasing in case 
(i) and strictly increasing in the other cases. In fact, for (i), (ii) and (iii) outer this 
follows already from (3.11) and [23]. As far as cases (iii) left and right are concerned, 
now the result follows from 

On the other hand, in all cases one has dk/difo ^ 0. Therefore, given T e I, the 
condition holds in order to determine ipo by the implicit function theorem so that the 
respective would have a period T. 

Finally, by using the above formulas one easily obtains that 

U v = -vd 2 x + v-2uj- 3{3(pl = va 2 [-d 2 y + 6k 2 sn 2 {y; k) + m] 

where y = ax and m = — 1 — 4fe 2 in cases (i) and (iii) outer, m = —4 — k 2 in cases 
(iii) left and right, m = — 1 — k 2 in case (ii), respectively. 

Lemma 3. The first five eigenvalues of the operator A defined by the differential 
expression A = — d 2 + 6k 2 sn 2 (y; k), with periodic boundary conditions on [0,4K(k)], 
are simple. These eigenvalues and their respective eigenf unctions are: 



fa 


= 2 + 2fe 2 -2Vl 


- fe 2 + fe 4 , 


My) 


= l-(l + fe 2 - y/l- 


- fe 2 + k A )sn 2 {y;k), 


/'I 


= 1 + k 2 , 




4>i(v) 


= cn(y; k)dn(y; k) = 


sn'(y; fe), 


1'2 


= l + 4fe 2 , 




fo(y) 


= sn(y; k)dn(y; k) = 


-cn'(y; fe), 


/'3 


= 4 + fe 2 , 






= sn(y; k)cn(y; k) = 


—k~ 2 dn'(y; fe), 


I'i 


= 2 + 2fe 2 + 2^1 


- fe 2 + fe 4 , 




= 1 - (1 + fe 2 + y/1 - 


-fe 2 + fe 4 )sn 2 (y;fe). 



Proof. The equalities Aip n (y) = faip n (y)i < n < 4 are established by calculation. 
By (5.5), (5.6) and the properties of elliptic functions, fa are simple and the rest of 
the spectrum lies in the interval (fa, oo). □ 

Corollary 1. The first three eigenvalues of the operator A, equipped with periodic 
boundary conditions on [0,2i^(fe)] ; are simple and equal to fa, fa, fa. 
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By Lemma 3 and its corollary, tp' v is the second eigenfunction of the operator 7i v in 
the cases: truncated pendulum and Duffing oscillator (left and right). For these cases, 
it makes sense to investigate the sign of d(v) which we do in the next proposition. 

Proposition 8. Assume that T <E I = (T , oo) and ip v = ip(x — vt), v > 0, is the 
periodic travelling-wave solution having a minimal period T . Then: 

(i) In the truncated pendulum case one has d{v) < 0. 

(ii) In the left {right) Duffing oscillator case, if 3v 2 — 8u 2 > 0, then d(v) > 0. If 
3v 2 — 8uj 2 < and 2v 2 — 2ujv — uo 2 > 0, then there is T max e (To, oo) depending only 
on the ratio uj/v, such that d{y) > for T > T max and d{v) < for T e (T ,T max ). 
If2v 2 - 2luv - cu 2 < 0, then d\v) < 0. 

Proof. We are going to use formula (4.9). In the example we deal with, one has 

2 TT( M 2k 2uJ ~ V 2 P 4 . 2U-V 2 3 4 
y =U {x,h) = — —X -—X, h=^—^ + -^. 

Given a nonnegative even integer n, denote I n {h) = § H=h x n ydx. It is well known 
that the linear space of integrals {I n (h),n even} forms a polynomial R[h] module 
with two generators, Io(h) and l2{h). Moreover, Iq and I2 satisfy a Picard-Fuchs 
system of dimension two. This implies that the ratio R{h) = I' 2 {h) / Ib(h) satisfies a 
Riccati equation. We shall use these facts to express d(v) as a quadratic form with 
respect to I' , I' 2 with polynomial coefficients in h and use the properties of the Riccati 
equation to determine the sign of d(v).The procedure might seem too long, but it is 
universal (at least in the case when a(u) and b(u) are polynomials) and therefore 
applicable to many other cases. 

Let us first express d(v) through the integrals I n . Below, we will denote for short 
the derivatives with respect to h by I' n) I" etc. Using the first equality in (3.13) with 
b = 0, we obtain 

x n dx 



H=h vy 



On the other hand, 

/ y 3 dx — (f U(x, h)ydx = '- — 1 { 

JH=h JH=h V 

By using these expressions, we obtain that 



2h T 2u-v 8 

h — 7T^- 

v 2v 




(5.7) 



For reader's convenience, below we proceed to derive the relations between integrals 
I n and the Picard-Fuchs system satisfied by Iq and I 2 . 

Lemma 4. (i) The following relations hold: 

{n + 6)pi n+3 + (2n + 6)(2uj-v)I n+1 = AnhI n - 1 , n = 1,3,5,.... (5.8) 
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(ii) The integrals I and I 2 satisfy the system 
AhI' Q - (2cu - v)V 2 = 3J , 



\r 0+ U + ^l)r^ 5h . (5 ' 9) 



3[3 u V 3/3 
Proof, (i). Integrating by parts, we obtain the identity 

/ [x n U'(x) + Inx^Uix^ydx = 0. (5.10) 

J H=h 

Indeed, 

/ x n U'(x)ydx = I x n ydy 2 = \l x n dy 3 = -hi I x n ~ x y 3 dx 

JH=h JH=h 6 JH=h 6 JH=h 

= -\n & x n ' 1 U(x)ydx. 

6 JH=h 

Identity (5.10) is equivalent to (5.8). 
(ii) Similarly, one has 

r x»U>(x)dx = r ^ / = _ 2n t ^ 

JH=h y JH=h y JH=h JH=h 

r x n U(x)dx r „ , 
& — — = & x n ydx. 

JH=h y JH=h 

Rewriting these identities by means of I n , we come to the formulas 



(5.11) 



(H'n+a + ( 2uJ - W ) J n+l = nI n-l, 

-(3I' n+4 - 2(2u - v)I' n+2 + 4hl' n = 2I n . 



(5.12) 



The last two relations imply 

Ahl' n - (2cu - v)I' n+2 = (n + 3)/ n . (5.13) 
Using (5.13) with n — 0, 2, we obtain the system 

Ahl' - (2cu - v)I' 2 = 3I , 
Ahl' 2 - (2uj - v)I' 4 = 5/ 2 . 

We remove I' 4 from the last equation by using the first equation in (5.12). As a result 
one obtains (5.9). □ 

After proving Lemma 4, we return to our main goal. Using the identities 

4/i 8(2lu-v) t 1 2m - v 
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and the first equation in (5.9), we remove J 4 , I' 4 and I from (5.7). The result is 



„ ( 3v 2 - Af3h , 2uj + 5v , 2u - v „ 



3f3v 2 



■ ir 2 J 2 



I" - -I'+L 



(5.14) 



Now, we differentiate (5.9) and determine the second derivatives from the obtained 
system: 

V(h)I% = — 4/i/q + (2uj - v)I' 2 , 



V{h)V{ = A(2UJ Q V \ l' Q + Ahl' 2 , 



(5.15) 



V(h) = 16h 2 + 



1 



h. 



Replacing in (5.14) and performing some direct calculations to simplify the result, we 
derive the formula 



d(v) 



1 



9vl 



8h 2 



12m 2 



h) Iq + (Au + I0v)hl' r 2 + (2v 2 - 2uov - uo 2 )I' 2 2 



By (5.15), the ratio R(h) = I 2 (h)/I' (h) satisfies the Riccati equation 
V(h)R' = A{2uJ ~ v \ + 8hR _ (2u - v)R 2 . 



(5.16) 



In order to remove the parameters from (5.16), we take new variables h, R through 



h 



8p 



(h + 1), R(h) 



2/3 



Thus we obtain the final formulas 



, (2u-v) 2 r 2 {h) - - 

d{v)= 72WW) W{k,R) 
where (below we will omit thoroughly the bars) 

w(h,R) = (4v 2 -4ujv-2L0 2 )R 2 + (4uj 2 + 8uJv-5v 2 )hR+(2uj-v) 2 h 2 + 3v 2 (R-h)-Quj 2 

(5.17) 

and R = R{h) satisfies the equation (and related system with respect to a dummy 
variable) 



4(1-^)7?' = l-2hR + R 2 , 



h = A(l-h 2 ), 
R=l-2hR + R 2 . 



(5.18) 



Since v > and Jo > (this is the area integral), the sign of d(v) is determined 
by w. The conic curve w(h, R) = divides the R 2 -plane into two parts W + and W- 
according to the sign of w(h, R). We have to identify the curve T in the phase portrait 
of (5.18) which corresponds to V 2 j I' and determine its location with respect to W + 
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and W_. In the case when d(v) changes sign, the main problem will be to prove that 
the conic curve intersects T at most once. 

It is easy to see that system (5.18) has two critical points (1, 1) and (—1, —1) 
which are saddle-nodes. They are connected by two separatrix trajectories, upper 
T u and lower Ti, please draw a picture. Obviously, the phase portrait of (5.18) is 
symmetric with respect to the origin. 

Lemma 5. In (5.18), the trajectory corresponding to the truncated pendulum case is 
Ft. The trajectory corresponding to the left and right Duffing oscillator cases is T u . 

Proof. If (fo £ A_, then ip^ G (0, V ~^ UJ ) in the truncated pendulum and right Duffing 

i 

(3 ' 0~ 



oscillator cases and tp\ G {ff^—jjr-, 2v g 4uJ ) in the left case. Correspondingly, one obtains 



h G ^0, — - — ^ - \ = (h c , h s ) in the truncated pendulum case, 

h G ^— - — ^ - , oj = (h c , h s ) in the left and right Duffing oscillator cases, 

where h c is the Hamiltonian level corresponding to a center and h s - to a saddle. 
Moreover, R(h) = I' 2 (h) / I'q(K) is analytic in a neighborhood of h = h c and, by the 
mean value theorem, 

lim R(h) = lim ^( = X 2 



^ Uh) x2 



We recall that X is the abscissa of the center and X = in the truncated pendulum 
case, Xq = v ~p U) in the left and right Duffing oscillator case. On its turn, by Picard- 
Lefschetz theory ([6], Chapter 3), near h = h s the integral I n (h) has the expansion 

I n (h) = I n {h s ) + a n (h - h s ) log \h - h s \ + (3 n (h - h s ) + 7„(/i - h s ) 2 log \h - h s \ + . . . , 

where «o 7^ 0. Therefore R(h) = I' 2 (h)/I' Q (h) is bounded near h = h s . In terms 
of the new coordinates (h,R), all this means that the phase trajectory of (5.18) 
corresponding to the truncated pendulum case is (h, R(h)) where —l<h<l, R{h) — > 
±1 as h — > ±1, and R(h) is analytic near h = —1. The unique phase curve with these 
properties is IV Similarly, the phase trajectory of (5.18) corresponding to the left 
and right Duffing oscillator cases has the same ends but is analytic near h — 1, hence 
it is T u . □ 

Below, we list some properties of the separatrices T/, T u and the conic curve 
w — 0. Apart of the separatrices, the conic curve depends on one parameter lj/v and 
undergoes several bifurcations when u/v varies. 

Lemma 6. The separatrices T u and Ti have the following properties: 

(i) T u is increasing and concave, Ti is increasing and convex. 

(ii) The tangent to T u at (—1, —1) is h = —1, the tangent to T/ at (1, 1) is h = 1. 

(iii) The tangent to T u at (1, 1) and the tangent to T/ at (—1, —1) have a slope \. 
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Proof. The proof easily follows from (5.18). Claim (iii) is a consequence of the ana- 
lyticity at the corresponding point. To establish (ii), we use the expansions containing 
logarithmic terms above. They imply that the asymptotic expansion near (—1,-1) 

of the non-analytic trajectories has the form —1 + /// \og(h + 1) + We determine 

pi from (5.18) and obtain that on T u one holds 

fl( ' ,) = - 1 -bi(lTT) +0 ( log " ( ' l + 1) )' (5 ' 19) 

A similar formula holds for Ti, hence (ii) is proved. Finally, in the strip \h\ < 1 one 
has 1 — 2hR + R 2 > 0, therefore R' > and all trajectories increase. Differentiating 
with respect to h the Riccati equation in (5.18), we determine the curve of inflection 
points R"(h) = in the phase portrait, which is (R-h)(R 2 + 2hR-3) = 0. As R' = \ 
on the line R = h, by (ii) and (iii) T u and Ti cannot intersect this line. Hence, T u 
lies in the domain corresponding to the concave trajectories, and Ti - to the convex 
ones. □ 

The proof of the statements listed below is obvious. 

Lemma 7. The conic curve w(h, R) = has the following properties: 

(i) It goes through the critical points (1,1) and ( — 1,-1) of (5.18) and w(h,h) = 
Qu 2 {h 2 -l). 

(ii) For u = it degenerates into v 2 (R — h)(AR — h + 3) = 0. For uj 2 + 2uv — 2v 2 = 
it degenerates into v 2 h[(R - 1 + (7 =F 4 v / 3)(/i - 1)] = 0. 

(iii) If(h,R(h)) are the local coordinates near (1, 1), then R'(l) = 1 - R"(l) = 

Suj i (uj 2 +2wv-2v 2 ) 
3d 6 

(iv) If (h(R), R) are the local coordinates near (—1, —1) and w ^ 0, then h\—l) = 0, 
h"{-\) = 2 " 2 - 3 2 J-" 2 . 

After the preparation done above, we proceed to prove statements (i) and (ii) of 
Proposition 8. To prove (i), let us denote by Q the open triangle in the (h, R) -plane 
having vertices at (1, 1), (—1, —1) and (1, — |). By Lemma 6, Ti C Q. On its turn, 
f2 C W 7 !. This holds because by Lemma 7 one has w(h, h) < for \h\ < 1, w(l, 1) = 
m(-l,-l) = 0, R'(l) < 1, h'(-l) = 0, and because w(l, -\) = ~u>(u> + 2v) < 
(recall that 2uj > v > in the truncated pendulum case). Proposition 8 (i) is 
established. 

The proof of (ii) is more complicated. In this case, 2u < v. Assume first that 
uj 2 + 2ujv — 2v 2 > 0. Then, by Lemma 7 (iii), (iv) and convexity, the conic curve is 
a hyperbola. One of its branches is contained in the half-plane h < — 1, the other 
branch lies above its tangent line at (1, 1), namely R — l — R'(l)(h — 1). As R'(l) < 
in the considered case, by Lemma 6 (i) and Lemma 7 (i) we conclude that T u C 
The same proof goes in the degenerate case u 2 + 2ujv — 2v 2 = 0. 

The case u 2 + 2uov — 2v 2 < is more delicate. If u — 0, then Lemma 6 (i) 
and Lemma 7 (ii) imply that T u C W + . Below we will take w / 0. Consider the 
branch R = r(h) of the conic curve going through both (1, 1) and (—1, —1). Writing 
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w(h, R) = 5 R 2 + 5iR + 5 2 , S > 0, one obtains 

r{h) = Z h±^R^ D = 5 2 1 -45 5 2 >0, \h\<l. (5.20) 
2d 

Below we will denote this branch by C. By (5.19) and Lemma 7 (iv) we have near 
h = -\ 



This yields that for ft, close to — 1, T u is placed above C. Similarly, by Lemma 6 (iii) 
and Lemma 7 (iii) one obtains that near h — 1, F u lies above C if 3t> 2 — 8cj 2 > 
and below C otherwise. Therefore, to finish the proof of Proposition 7, we have to 
establish the following: 1) T u is entirely placed above C if 3v 2 — 8uj 2 > 0, and 2) 
r u intersects the conic just once if 3v 2 — 8uj 2 < 0. In the first case, we would have 
T u C W + , while in the second one, the part of T u near h = 1 would be in W_ and 
the remaining part - in W + . Unfortunately, both T u and C are concave and it is not 
so easy to determine the number of their intersections. 

Below, we proceed to determine the number of contact points that C has with the 
vector field (5.18). Because of the type of critical points, in our case the number of 
intersections is less than or equal to the number of contact points. As well known, 
the equation of the contact points is given by 

j- g {R- r{h))\ R=r(h) = [R- hr'(h)] lR=r(h) = 1 - 2hr(h) + r\h) - 4(1 - h 2 )r\h) = 0. 

(5.21) 

We replace r' = —(5[r + 5' 2 )/(25 r + 5i) in (5.21) and use once again the quadratic 
equation satisfied by r to obtain 



8 1 + A5 2 h+ 5 -4 1 + 4^(1 -h 2 ) 



_l_ r 



2 



s 2 

25 + 2<?i/i - 25 2 + -± + 461(1 - h 2 ) 



0. 



Next, replacing r from (5.20) and performing the needed calculations we get 
[125 G {2uv - v 2 )(l - h 2 ) + D}VD - [3v 2 (l + h)D + 25 (1 - h 2 )D'} = 0. 
By (5.20) and (5.17) one has 

D= (l + ft) 2 (9t; 4 + 24^ 2 5oC), D' = (1 + h)(18v A + 24lu 2 5 (( - 1)), C = \— \ > 0. 

We replace these values in the equation above and divide the result by (1 + h) 3 . One 
obtains 

[9w 4 + 12(2u 2 + 2uv - v 2 )5 (W9v 4 + 24u}' 2 5 ( 
-[27v 6 + 12{3v 4 + 6u 2 v 2 - 4u 2 5 )5 o ( + 48uj 2 5 2 ( 2 } = 0. 
At the end, we introduce a new variable z by 



^9v 4 + 24uj 2 5 ( = 3(v 2 + z), z £ (0, oo) 
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to obtain the equation of contact points in a final form P(z) = where 

P(z) = 3z 3 + 6(3v 2 - 2ujv - 2u 2 )z 2 

+A(v + cu)(9v 3 - ISujv 2 + Aoj 2 v + Auj 3 )z + 2v 2 5 (3v 2 - 8cu 2 ). 

Taking P(z) = 3z 3 + A 2 z 2 + A\z + A , it is easy to verify that: 

A 2 >0,A 1 >0,A >0 for 3v 2 - 8cu 2 > 0; 

A 2 > 0, Ai > 0, A < for 3v 2 - 8cj 2 < 0, v + tu > 0; 

A 1 < 0, A < for v + to < 0; 
(recall that f — 2uj > and 5 = Av 2 — Auov — 2oo 2 > 0). By Descartes chain rule, 
the equation P{z) = has no positive root if 3v 2 — 8u 2 > and has exactly one 
positive root if 3v 2 — 8u 2 < and 2v 2 — 2uv — u 2 > 0. All changes of the variables 
we used throughout the proof were one-to-one, hence there is no contact point in the 
first case and there is just one contact point in the second case. As a result, in the 
first case C does not intersect T u and therefore T u C W + . In the second case, C 
intersects Y u at a unique point corresponding to some h G (—1,1) so that the part of 
T u related to (—1, ho) is in W + and the remaining part is in W-. Let y?o be the value 
corresponding to ho according to (3.11) and let T max be the period of the orbit going 
through (</?o, 0) in the (</?, </?')-plane. Then the orbits from the period annulus having 
a period T > T max will be stable and the remaining ones unstable. Proposition 8 is 
completely proved. □. 

Let us recall again that in the left(right) D uffing oscillator case the period T belongs 
to the interval / = (T , oo) where T = 27ry / ^3^. Lemma 3 and Proposition 8 imply: 

Corollary 2. Both conditions of Assumption 1 are satisfied in the left(right) Duffing 
oscillator case, provided that: 

(i) 3v 2 - 8tu 2 > 0; 

(ii) 3v 2 — 8uj 2 < 0, 2v 2 — 2ojv — u 2 > and the period T is sufficiently large. 

Proof of Theorem I. Theorem I is a direct consequence of Corollary 2 and Theorem 
2, taking into account that <p does not oscillate around zero only in the left and right 
Duffing oscillator cases. □ 

Example 5.3: Coherent single power nonlinearities. Let a(u) = (3u k+1 , b(u) = 
7/9« fc , k e N, where |7| < 1 and v/f3 > 0. Or more generally, one can take a(u) = 
2uu + j3u k+1 , b(u) = 2uj + / ~f/3u k , where r = (v — 2u)//3 > 0. For simplicity, assume 
that c\ = 0. Then (3.5) is satisfied with X = r 1 ^ and Proposition 1 works. Next, 
equation (3.3) implies that ip depends on r alone. Moreover, one has 

<p = r 1/k ^, <p = r 1/k (p , </?i = r 1/fe £>i, c 2 = r^c 2 , U(<p) = r*U(<p), 
where <p, <^o; <^i, c 2 and U do not depend on v, uj and j3. Therefore we obtain 
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= w^T) L <^ 2 + ^ = w^d /. + ^ )<ic - (5 - 22) 

Therefore d(f ) takes the sign of (3 (and v — 2uj). 
6. The perturbed gBBM equation. 

Let us denote by C(X, Y) (respectively, B(X,Y)) the set of all closed (respectively, 
bounded) operators S : X — > K. When X = F, we shall write simply C(X) and £>(X). 
USE C(X,Y), denote by G(S) C X x F its graph. Below we choose X = Y = 
L 2 [0, T}). One can define a metric 5 on C(L 2 [0, T]) as follows: for R, S E C(L 2 [0, T]), 
set 

8(R,S) = \\Pr — Ps\\b(l 2 xl 2 ) 

where Pr and Ps are the orthogonal projections on the graphs of G(R) and G(S), 
respectively. 

The following results are well known (see Kato [M], Chapter 4, Theorems 2.14 
and 2.17): 

Theorem A. Take X = Y = L 2 [0,T] and assume that R, S G C(X,Y), A 6 B(X,Y). 
Then 

5(R + A,S + A) < 2(1 + \\A\\ 2 )5(R,S). 

Theorem B. Assume that S G C(X,Y) and B is a S-bounded operator satisfying 
\\Bf\\ < a||/|| + 6||5 , /|| with b<l. Then 

R = S + B G C(X, Y) and 5(R, S) < 

l — o 

Let us introduce a small perturbation in (gBBM) by taking a{u) = 2uu+j3u k , k = 2, 3 
and b(u) = ^g{u) in the general equation (1.1) where 7 is a small real parameter. We 
will denote the travelling-wave solution corresponding to (1.1) again by (p v and the 
one corresponding to (gBBM) (that is, when 7 = 0) by ip° v . If T is the period ip®, 
then by Proposition 3 for sufficiently small 7 there is a smooth function ^0(7) such 
that the wave solution ip v of the perturbed equation which goes through the point 
((^0(7), 0) in the (<p, <//)-plane will have the same period T. For the related operators 
Ti v we have respectively 

H v = -vd 2 x +v- a'{ip v ) + i[g(vv)d 2 x + g'M<p' v d x + \gf\<p v )<p% + g'&vWX 
H° v = -vd 2 x +v-a'(ip° v ). 

Theorem 3. If Assumption l(i) holds for the operator Ti®, then for 7 sufficiently 
small it also holds for the operator 7i v . 

Proof. In order to apply Theorem A, we define the operator A to be a multiplication 
by the function —a'((p®) and take S = —vd 2 + v. Then denoting £> 7 = H v — TC®, we 
obtain from Theorem A the estimate 

6(H V , H° v ) = S(S + + A, S + A) < 2(1 + | \A\ \ 2 )5(S + S 7 , S). (6.1) 
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As 

B 7 = a\ifl) - a'{<p v ) + 7 [-tr 1 ( ? (y^)S + g'M<p' v d x + G] 

where 

G = g(<Pv) + ^"{VvWv + 9'M<p", 
for / G D(S) = Hp er , we further have 

||5 7 /|| <m^\a'(^)-a'M\.\\f\\ 

+ \l\[max\g(<p v )/v\.\\Sf\\ + max\g'{<p v )<p' v \.\\d x f\\ + max|G'| 

From Plancherel identity, we have 



\Sf\\ 2 = [ Sf.Sfdx = T £ \Sf{k)\ 

JO L._ „ 



2 



k=—oc 

+00 +00 

= v 2 E (l + \k\ 2 ) 2 \J(k)\ 2 >K(v) £ \kmk)\ 2 = K(v)\\d x f\\ 2 

k=—oo k=—oo 

where K(v ) is a constant depending only on v. By using Propositions 3 and 4, we 
obtain that max |a'(y9°) — a'((p v )\ = 0(7). Taking into account these estimates, we 

come to the conclusion that the operator B = B 1 is S"-bounded, with constants a, b 
tending to zero as 7 — > 0. Applying Theorem B and inequality (6.1), we obtain 

5(H V ,H° V ) < 2(1 + ||A|| 2 )^±? = 0( 7 ). (6.2) 

Therefore by [31] (Theorem 3.16 in Chapter 4), the operators 7i v and 7i° v have the 
same spectral properties. More precisely (see also |Tj, p. 363), let U is the open disc 
bounded by Y and speciH®) DU — {a ,0}, where a is the negative eigenvalue of 
7i®. Choose circular contours T\ and T2 contained in U, such that if U\ and U2 are 
the open discs bounded by Ti and T 2 respectively, then spec(?i®) D U\ = {«o} and 
spec(7^°)nf/ 2 = {0}. From (6.2) for sufficiently small 7 we have that the spec(H v )r\U 1 
and spec(TC v ) C\U 2 consist of a single, simple eigenvalue. Since is an eigenvalue of 
Tt v we must have spec(H v ) fl U 2 = {0}, which shows that is a simple eigenvalue of 
TL V . Similarly we obtain that spec(H v ) fl U consists of a finite set of eigenvalues of 
total multiplicity 2. This completes the proof of the theorem. 

Proof of Theorems II and III. Take u — 0. Inequality (5.22) implies that d(v) 
is positive. When 7 = 0, applying respectively Proposition 7 and Corollary 1, we 
conclude that TL® has the needed spectral properties in both cases. By Theorem 3 
the same remains true for small I7I as well . Hence, Theorem 2 applies to both cases. 
□ 

7. Small-amplitude travelling-wave solutions. 
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Since the expressions in (3.10) and (4.9) are too complex, it is difficult to use them 
in general. Below we are going to consider the simpler case when the periodic waves 
we study are of small amplitude. That is, ipi — ip is close to zero and therefore 
the periodic trajectory is entirely contained in a small neighborhood of some center 
(X , 0) G M 2 given by (3.5). Below we establish that an important role in the stability 
of the small-amplitude travelling-wave solutions is played by the first isochronous (or 
period) constant. Let us recall that the period function has an expansion 

T(r) = T + T 2 r 2 + T 4 r 4 + T 6 r 6 + T 8 r 8 + . . . 

with respect to r, the distance between the center at (X , 0) and the intersection point 
of the orbit with the x-axis ((p ,0). Then T 2 k, k = 1,2, .. . is the kth isochronous 
constant. When all period constants are zero, the center is isochronous and all orbits 
surrounding it have the same period T . Since T'({p ) = —dT/dr, the above expansion 
yields immediately: 

Proposition 9. The period of a small- amplitude travelling-wave solution around a 
non-isochronous center is not critical. 

Let us mention, however, that the calculation of the isochronous constants is a difficult 
task and the problem of isochronicity is completely solved for only few particular cases. 
Denote for short 

a"(X ) _ a"'(X ) b'(X ) b"(X ) 

P 3(a'(X )-vY q 3(a'(X )-vY m v - b(X )' " v - b(X ) ' l ' J 

As we shall see below, the first isochronous constant in our case is expressed by 
T 2 = T tf, T = 2tt 



v - b{X ) M _ 1 2 



\ l(v v , ^ = T^( 15 P +3g-m -Qpm-2n). 

\ a'(X ) — v 16 

Take a small positive e and let (p = X — e. Then by (3.9) we obtain 

<P! =X + e + pe 2 +p 2 e 3 + O(e 4 ). (7.2) 

The expression (7.2) as well as the formulas that follow are obtained by long and 
boring asymptotic calculations which we will omit here. Thus, by using (3.8) we 
come to the expression 



U(X) 



eV(X ) -y)(l- s 2 )[l + p(l - sje - \{p\l + 6s) + qjl + s 2 ))e 2 + ..] 



(v - 6(X ))[1 - mse - \(mp(l + s) + ns 2 )e 2 + . . .] 

Then, by (3.7) and (3.8) we obtain 

T = T [l + ^ 2 + O(e 4 )]. (7.3) 

In order to calculate the sign of the derivative d(v), we prefer to apply a more direct 
approach, deriving first an alternative formula instead of (4.9). By using (4.3), we 
obtain 
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fo yU(s 



In the calculation below, we denote by ip , ipi etc. the derivatives with respect to the 
parameter v. As U(s) vanishes at ip and ipi, then clearly 

d ( m r— \ m (U(s) + ipoU^s^ds 



i(i:^)=L 



/ J 90 2yJU{s) 

To perform differentiation in the other part of the integral, we first use a change of the 
variable like (3.8) and return to the initial variable after the calculations are made. 
Thus, 



d ( fPi s 2 ds \ (fi — (po f^ 1 s 2 ds 
dv 



+ 



• sju{s)) Vi-WJvo ^U(s) 
vi (0i - (pp)s + v? y?i - y?iV?o 2sds 
w <Pi ~ V?o yJU(s) 

f 1 (<Pi - 0o)s + <Po<Pi -<Pi<Po TTl , ^\ s2 ds 



f Vl (tt( \ j ■ tt n i (0i -yo)s + Wi-^o rr ,, " 
-/ [U(s) + (p U V0 {s) + U(s) 



As a result, one obtains 



<i(t,) = l wkkry {7A) 



with 



V( 8 ) = 2 ^^s 2 U(s) + (U(s) - s 2 )(U(s) + tpoUM) 
tp 1 -(p 



+ (AsU(s) - s 2 U\s)) 



2tt'/ \\ (V^l - 0O)S + VWl - 01^0 



(7.5) 



<Pi - V?o 

From (3.5), (7.1) and (7.2), for the derivatives with respect to v we get respectively 

= Xo ° (p + 3p j + x ^ + 0^ + + ^ + + (7 ' 6) 

Replacing (7.2) and (7.6) in the expression of V(X) given by (7.5), one obtains 
, (l-s 2 )(a'(X )-v)e 2 e 



v - b(X ) 



( P X 2 + 4X -mX 2 )s 



+ + ^ + Xoip(4 + 8s- 2s 2 ) + 2ms 2 } 

\ v-b(X ) 

+X 2 [p 2 (3 + 2s) + pm(2s 2 - 2s - 1) + |(1 + s 2 ) - (2m 2 + n)s 2 ]^j e + 0(e 2 ) 
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The leading term of e can be determined from (7.3) after differentiation with respect 
to v. One obtains 



e ~ 



2Tk " 4(a'(X ) - t;)ife r ° m " 3PX ° " 1 " v - b(X ) , 



Finally, replacing in the integral (7.4), we derive the formula 



d(v) = (a'(X ) - v) 



t 2 



(a'(X ) - v) 2 2T 2 



+ 0(e). 



(7.7) 



Example 1: The BBM equation. Let a(u) = 2uju + |u , 6(m) = 0. Assume for 
definiteness that t> > 0. The other case is considered similarly. Using the notation 
introduced earlier and the expression of X , we obtain 



i\X )-v = D 1 '\ p = 



1 15 

g = m = n = , ^ = y^, T = 2vr 



£1/4 > 



To 



15T 
16L> ; 



T (2o; 2 - aw + 3ci) 

^73 



, d( v ) 



DV2 



2 8(2o; 2 - uv + 3c^ 2 



15v 2 



Now, expressing ci and X through D, we rewrite the expression in the brackets of 
d(v) as 



A 



, 2(D - v 2 + 2ujv) 2 



15v 2 



15v 2 

x 



D 



D 



vD 1 ' 2 



)v(v - 2u) 



X 



vD i/2 _ (1 + _ 2w ) 



Next, denote by Ai, A 2 the roots (with respect to D 1 / 2 ) of the first multiplier, and 
respectively by A3, A4 the roots of the second multiplier. 

a) If < v < ^-(9 + 5^/|)a;, then both multipliers are positive and A < 0. 

b) If ^-(9 + 5^/|)u; < v < 2a;, then the first multiplier is positive, < A4 < A3 and 
therefore A > <^> A 4 < D 1/2 < A 3 . 

c) If A > 2u, then < A : < A 3 , A 2 < 0, A 4 < 0. Therefore A > <£> Ax < D 1 / 2 < A 3 . 
Taking a second power of the above inequalities, it follows that 



d{v) > <£> 



\ 2 

A 4 



2a;) 2 < 6ci < A 



2a;) 2 in case b), 



A 2 — (v — 2uj) 2 < 6ci < A 2 — (v — 2a;) 2 in case c). 



(7.8) 



As Afe are homogeneous of first degree with respect to v,u, one can reformulate (7.8) 
as: d(v) > if and only if (cu/v,Ci/v 2 ) £ Q C 1R 2 where Q can be explicitly written 
down if needed by using the formulas of the quadratic roots A&. We are not going to 
do this. 
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Example 2: The modified BBM equation. Let a(u) = 2uu + (3u 3 , b(u) = 0, 
ci = 0. 

(i) X = (Global center or truncated pendulum case). Then d(v) = -j^Tq + 0(e) 

(ii) Xl = (Duffing oscillator). Then d(v) = %=^T +O(e) and = > 0. 

As seen from these examples, d(v) could be negative or positive, depending on the 
case. 

Proof of Theorems IV and V. In Examples 1 and 2 above, we calculated the 
first term of d(v) provided that b(u) = and understood when it is positive. If one 
take b(u) = 7g(w) with |7| small, then d(v) will differ from the case b(u) = by 
a term which is 0(7). Therefore d(v) will keep the sign of its first term provided 
that I7I and the amplitude e are small enough. This means that d(v) > in the 
domain Q from the first example and for 3v 2 — 8u 2 > in the Duffing oscillator case 
of the second example (recall that all solutions not oscillating around zero belong 
to this case). Since the operator related to the above cases satisfies the spectral 
properties required in Assumption 1, the statements of Theorems IV and V follow 
from Theorems 2 and 3. □ 
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